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Abstract 

(-) 
From a time-dependent analysis of the decay B°—> ir + ir and using a 

model-dependent A P /A T ratio of the penguin-to-tree amplitudes contri- 
buting to the decay, both the weak phase a and the strong phase shift 
difference 5 can be extracted from the data. The value of the weak phase 
P, expected to be measured from the decay B° — ► J/ipKg, is used to pa- 
rameterize the value of A P /A T and the corresponding penguin correction 
to the observed asymmetry in B° d — > tt + it~ . 



1 Introduction 

(-) , _ (-) 

Measurements of CP violation in -B°— > n + n and B° d -+ J /ipK° s decays are ex- 
pected to provide information, respectively, on the angles a and /3 of the unitarity 
triangle [l|] . As it is well known, in B° — > J/ipK° s the weak phase of the penguin 
contribution is the same as in the (largely) dominant tree amplitude and there- 
fore the value of sin 2/3 can be directly extracted from the measurement of the 
time-dependent CP asymmetry with no hadronic uncertainties. In particular, as 
being dominated by a single amplitude, this channel should be insensitive to the 
presence of phases of strong-interaction origin. 

On the other hand, the extraction of sin 2a from a time-dependent analysis 
of the asymmetry in B° — > it + it~ appears to be less straightforward 0. This 
is mainly due to the presence of gluonic and electroweak penguin terms with 
weak phases different from those of the dominant tree b — > udu amplitude. The 
common prejudice is that the penguin should be smaller than the tree diagram. 
However, for a quantitatively meaningful statement, the relative size must be 
either estimated theoretically or measured experimentally from a suitable data 
sample of -B-decays. An additional source of uncertainty is introduced by strong- 
interaction phases which can be different for tree and penguin amplitudes and 
difficult to model at the non-perturbative hadronic level. 

While the size of the tree and penguin amplitudes (in particular, their ra- 
tio) can be somehow assessed under some theoretical assumptions, e.g., the fa- 
ctorization hypothesis for non-leptonic decay amplitudes 0, the effect of the 
strong-interaction phases is much harder to estimate reliably ||. The conven- 
tional expectation is that such phases should not be large. This seems verified 
by the phases of perturbative origin 0,0, coming from the cc penguin with 
on-shell cc quarks |7|, M, and is expected to be true also for the final state, long 
distance, non-perturbative phases, considering the high velocity of the produced 
7r + 7r~ pair. Nevertheless, from considerations based on Regge phenomenology, it 
has been recently suggested that, contrary to expectations, final state interactions 
might be important even in B non leptonic decays |§. Also, significant effects 
have been estimated by semiquantitative (and model-dependent) calculations for 
the color suppressed channels B° — > D°7i° and B° — > 7r°7r° JH], |l ij] . 

Several strategies have been proposed to extract the value of sin 2a in a model- 
independent way avoiding the difficulty associated with penguins in B° — > 7t + tt~ . 
However, although quite effective in principle, their application to a practical 
analysis of the data might run into some difficulties and the effective potential 
of these methods have still to be realistically assessed by taking into account the 
expected statistics and experimental sensitivities. For example, the isospin anal- 
ysis of B° — > 7T7T decays |L2|, |13|],[[L4]] which would enable to directly reconstruct 
the unitarity triangle, suffers from the predicted low rate for B° d — ► 7r°7T , of the 
order of 1CT 6 or less fl5| , 16|, and unfavourable background conditions. Alter- 
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native methods require more complicated analyses to fit penguin amplitudes by 
combining measurements of B Q d — > 7t + tt~ with different penguin-dominated pro- 
cesses such as B -»• Ktx O, M, III EO, ETD, E3 or £° -»■ # i?° E3] and relying on 



SU(3) symmetry of matrix elements and/or first-order SU(3) breaking. In that 
case, the sensitivity to a could be rather difficult to fully establish in practice. 

The hadronic uncertainty in the extraction of sin 2a from the CP asymmetry 
(-) 
in B°-^ 7r + 7r has been discussed by many authors and the analysis of this channel 

extended to other, non-strange, AS* = final states such as pir or pp, and the 

corresponding AS = 1 decays [|5], |23], |25|, ^B|, E7J. To derive quantitative results, 



assumptions are made on the relevant hadronic matrix elements (in general, the 
factorization model), and on the final state strong interaction effects which in 
most cases are assumed to be negligible . 

In the following, we focus on the determination of sin 2a from the measured 
(-) 
CP time asymmetry in B°—> n + n . A discussion of the penguin-induced correc- 
tion Aa in terms of a itself and of the strong phase 5 = 5 T — 5 P between tree 
and penguin amplitudes was presented in @, where a first-order expansion of the 
asymmetry in the penguin-to-tree amplitude ratio A P /A T was used. As a result 
of this analysis, it turns out that Aa might be significant, except for large values 
of the asymmetry and that, in any case, the correction should be maximal for 

s = o. Q 

The strong phase effect on Aa for general values of S has been further con- 
sidered in H using the exact expressions for the CP asymmetry and the tree-to- 
penguin amplitude ratio from the factorization model as input. The indication 
from |J is that, within the assumed model, the expected relative shift in a should 
not exceed 30% and, similarly to 0, is maximal for 5 = 0. 

In this note, we take a somewhat different point of view and try to gain in- 
sight into the actual value of the strong phases difference 5, without a priori 
assumptions on their size, assuming that both the mixing-induced and the direct 
CP violation terms will be measured from the time-dependent CP asymmetry of 

B°— ► 7r + 7r . Similarly to most previous analyses, we assume the approximation 
where the top mediated penguin dominates (in which case the electroweak phase 
of the penguin can be identified with the angle (3) and adopt the model dependent 
estimate of the ratio A P /A T provided by factorization. This ratio then depends 
on both the "true" value of a and on (3. Using this relation, we evaluate the cor- 
rection to the "measured" value of a (obtained from the fit of the time-dependent 

(-) 
rates in B\—* n + n ), required to determine the "true" value of a, as a function 

of /3 alone. For this correction we use the exact expression for the penguin con- 
tribution to the time asymmetry, so that the analysis potentially applies also to 
the case of large penguin-to-tree amplitude ratio. We also gain information on 



1 Related analyses can be found, e.g., in |2£ 



the effect of the strong phase S, corresponding to the different possible values of 
(3. To summarize, from a time-dependent analysis of B° d — + 7r + 7r~, combined with 
a clean and accurate independent measurement of the angle j3, we show that one 
can in principle: 

• extract the measured uncorrected value of a 

• evaluate (a model dependent) correction Aa 

(-) 

• gain insight into the role played by S in B° d — > 7r + 7r . 

The price is the introduction of a modest, model-dependent, correlation between 
two otherwise independent measurements, and the estimate of the ratio A P /A T 
from factorization. 



2 Model-dependent evaluation of A P /A T 

To make the presentation self-contained, this section briefly reviews the basic 

(-) 
points of the derivation of the penguin-to-tree ratio for B°— > n + n in the facto- 
rization approach. 

The relevant effective weak Hamiltonian can be parameterized as [^, |30], [31 

AG 



H, 



eff 
W 



V2 



v ub v: d (ciOi + c 2 o 2 ) - v tb v; d £ aOi 



+ h.c, (1) 



where c, are short-distance Wilson coefficients defined at a scale \x of the order 

of the heavy quark mass m b , and O; are a set of local quark operators with the 

appropriate quantum numbers. The first two terms in ([I]) represent the tree 

diagrams while the other ones are the contributions of strong and electroweak 

penguin diagrams. For simplicity we can neglect electroweak penguins, which 

(-) 
are found to give a small contribution to the mode -B°— > n + n of interest here 

H| [31], |32|]. In any case they could be included following the remarks in |[26|| . 



Retaining only strong penguin operators, the sum in (|T|) reduces to i < 6 and the 
explicit expressions of the relevant O, are: 

Ox = d a j^Lu p u^^Lb a ; 2 = dj^Luu^Lb; 

3 = d^Lb E q > qYLq'; 4 = d a ^Lb? E q > q'^Lq'*; (2) 

5 = d 7ll Lb E q > qYRq'; 6 = d^^Lb? ^ q'^Rq' a , 

where a, (3 are color indices; L, R = | (1 =)= 75); and q' runs over all quark flavors. 
For the corresponding Wilson coefficients at the scale fi ~ rrib = 4.8 GeV we use 
the values [BlU 



ci = -0.315; c 2 = 1.150; c 3 = 0.0174; , . 

C4 = -0.0373; c 5 = 0.0104; c 6 = -0.0459. ( ' 



In the factorization hypothesis, after Fierz reordering the operator Oq, one 
directly obtains from ([!]) and (0) 



(K + K-\H%>\B°) = i^ 



V2 



m B — mz 



I) UF^{ml)]{v ub V: d (jjC 1 + c 2 



-v tb v t * d 



(j?c 3 + c 4 ) + (^c 5 + 



C6 



2lllZ 



} 



(4) 



(m b -m u )(m d +m u ) _ 

where N is the number of colors. The following definitions have been used: 

_ _ i 

(7r~|d7 M L«|0) = -(7r~\djf,Ru\0) = -f n (p n ) 



(5) 



with /„- = 132 MeV the pion decay constant, so that from quarks equations of 
motion 

i rnlU 



(7i-\dLu\0) = -(ir~\dRu\0) = -- 

2m d + m u ' 



(6) 



and for the B° — > n + l vi matrix elements |]33 



(7i+\u lfl Lb\B°) = (n + \u lfl Rb\B°) 



(Pb + Pn)^ + 



'-% 



pB-^tf) + \^^ %F B^ {(f) 



(7) 



where q = ps — Pn, so that from equations of motion 
(<K + \uLb\B°) = (7i + \uRb\B°) = - 



\m\-m\ B _. 7i ., 



2 m b 



m„ 



{<?)■ 



(8) 



With both the coefficients q and the matrix elements of Oj real, the penguin 
amplitude in Eq. (^) has the (weak) phase of VtbV t * d . It is well-known that, 
by the unitarity of the CKM matrix, the top penguin dominance occurs if the 
difference between the c- and li-mediated penguins can be neglected, which is the 
case for rn^/ml ~ 0. However, phenomenological estimates indicate the values 
m c = 1.3 — 1.5GeV and rat, = 4.5 — 5GeV [0, |35[] . Thus, the verification of 
such approximation is clearly a point deserving further analysis.^ As a matter 
of fact, the consistent treatment of both the Wilson coefficients and the matrix 
elements of operators at next-to-leading order in QCD provides modifications of 
order a s to the simple structure of Eqs. (Q) and (|]). In particular it introduces 
imaginary parts into the coefficients [|T^, |3T| which determine final state strong 
interactions phases at the perturbative quark level. The actual value of such 
phase is somewhat parameter-dependent, although generally quite small, so that 



2 An attempt of extracting sin 2a from B° d — > 7r + 7r~ relaxing this approximation has been 
recently discussed in |22j ], assuming the approximate linear expansion in A P /A T of the CP 
asymmetry and the SU(3) relation to the B — > itK rate. 



in the following we will consider the strong interaction phase as a free parameter 
with a priori arbitrary possible values including (but in principle not coinciding 
with) the perturbative effects. 



(-) 



The general decomposition of the decay amplitude for B° d — ► 7i + n in terms of 
tree plus penguin contributions reads: 



A 
A 



A(B° 
A(B° 



7T+7T-) = A T e i{5T+M + A P e i{5p+ ^ p) 
7T+7T-) = A T e i{&T -^ T) + A P e i{ - 8p -' t,p) , 



(9) 



where (0r, 4>p) an d ($t, Sp) are weak and strong phases, respectively. In the top 
dominance approximation, (0) can be rewritten as 



A = \V: b V ud \e^ T+ ^T+\V: b V td \e^ p -^P 
A = \V ub V: d \e^ T -^T+\V tb V td \e^ 5p+ ^P, 



(10) 



where the CKM entries have been factored from the pure hadronic matrix ele- 
ments T and P. In the leading order Wolfenstein parameterization, the unitarity 
triangle angles : 



a = arg 



v td v t t 



(3 = arg 



\ VudV* b J \ 

are related to the parameters p and rj as 

V 






7 = arg 



v cd v; b , 



tana 



7] 2 - p(l- p) 



tan/3 



tan 7 



1 
P 



;n) 



(12) 



Comparing ( [TOD to (|j) one directly reads the expressions of T and P in the 
factorization approach: 



T 



V2 



(ml - ml) UFf+*(ml) 



N 



c\ + c 2 



(13) 



P 



x 



Grp 

V2 



m B -mAUF 



lS ^(m!' 



^ 3 + Q ) + (^ C5 



Q> 



2m? 



(m b - m u ) (m d + m u ) 



(14) 



From ( [T3"D and (|14[) one can notice the property of the ratio P/T of being in- 
dependent of particular models for the B — > ir form factors. Thus, apart from 
the factorization hypothesis itself, it only depends on the short- distance coef- 
ficients which can slightly vary according to the choice of the scale /i, and on 
the current quark masses in the last term of flUp. We choose m b = 4.8 GeV, 



m-d + fnu = !5MeV so that, with the coefficients ([3]) and N = 3 (the dependence 
on N is weak), one finds, as in || |26|, P/T ~ 0.055 . Reasonable variations 
of quark masses around the chosen values, taking into account the spread of 
the various determinations, can affect the term mentioned above by about 30%. 
However, being anyway suppressed in ( ^4|) by the large rrib and small short dis- 
tance penguin coefficients, this reflects into an uncertainty on the ratio P/T of 
the order of 10%. 



3 Discussion on the penguin correction Aa 

The time-dependent decay rates are given by : 

I "J 12 _|_ I A 12 

T(B°(B°) -> 7T+7T-) = l - l —-^ L e - r l*l (1 ± a sin Ami ± b cosAmt) (15) 

where the two real coefficients ao and bo are : 

21m(e 2if3 AT) 
a = 



| A | 2 + | A | 2 

(16) 

A I 2 - I A I 2 



A| 2 + M 



In terms of the ratio of the penguin to tree contribution A P /A T , the two coeffi- 
cients a and b$ can be written as : 

— sin la + 2 (A P /A T ) cos 5 sin a . . 

a ° '' = 1 - 2 (i4 P /A T ) cos S cos a + {A P /A T ) 2 ^ ' 



bo 



—2 (A P /A T ) sin 5 sin a 



1-2 (A P /A T ) cos <5 cos a + (A P /A T ) 2 



where 5 = 5 T — 5 P is the strong phases difference and, as in eq.([H]), T — (fi P = 
7 + /3 = 7r — a in the SM. When A P /A T = , the b term vanishes and a is 
extracted from the simple relation ao = — sin 2a . Else, a Aa correction term 
which depends both on a and cos 5 , has to be added. The Aa correction is 
therefore maximal for 5 = 0, as pointed out previously. 

The ratio of the CKM matrix elements in A P /A T can be expressed in terms of 
the angles a and (3 of the unitary triangle as : 

I MS I = ™1 (19) 

\V ub V: d \ sin/3 y ' 



and therefore, in terms of the ratio of pure hadronic matrix elements P/T 

P. 



J\.rp 



(sin a cot j3 + cos a) 



T' 



(20) 



The expected value of the ratio A P /A T — calculated using the P/T value obtained 
in the previous section — is plotted in fig. [J] as a function of a for four different 
values of /3. The curves in fig. [TJ are drawn for /3 = 11°, 18°, 25°, 32°, respectively 
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Figure 1: A P /A T plotted against the weak phase a for four different values of 
(3 = 11° (dashed upper curve), 18° (solid), 25° (dashed), 32° (dotted) . 

and the highest (lowest) curve lies close to the present experimental lower (upper) 
limit on (3 |]3"f5|],|]3"7|],|3"8|]. In the same picture, the current experimental bounds on 



a are superimposed onto the horizontal axis. 

From inspection of fig. p], we would conclude that the value of A P /A T predicted 
by the model could be as large as 0.3 if the value of f3 would turn out to be close 
to its present lower limit. However, we must keep in mind that the allowed ranges 
for sin 2a and sin 2(3 are correlated through the common dependence on p and rj 
(see figs. |2](a),(b)). Therefore, the values of a and (3 are allowed to cover a subset 
of the region bounded in fig. p] by the envelope of the (3 curves and by the a limits 
on the ascissa. This is best shown in fig. 0(c) where the present limits on the 
CKM triangle |3^|,[^7j,||38| are used to bound the shaded area which represents 



the allowed value of A P /A T as a function of a. From the picture, we see that the 
penguin to tree amplitude ratio can, at present, vary from a minimum of 0.08 to 
a maximum of 0.26, approximately. 
From equation (|17]), the correction Aa to the observed value of the a^ coefficient 

Aao = clq + sin 2a (21) 

is plotted in fig. |||(a) for 5 = 0, as a function of a. For the same value of 5, 
fig. |3](b) shows the Aa correction obtained, as in 0, by expanding equation 
( |17D to first order in A P /A T 

Aa = —2 (A P /A T ) sin a cos 2a cos 6 (22) 

and taking 5 = 0, while fig. §(c) shows the difference between the exact and the 
first order result. A common feature is that the Aao correction is positive in 
the interval tt/A < a < 37r/4 while it has the opposite sign and is smaller in 
magnitude, outside this range. However, while the first order expansion predicts 
a correction which vanishes for a = ir/A and a = 37r/4 ( for any value of (3 ) and 
is maximal for a = n/2, the second order term introduces a dependence on /3 in 
the position of the zeroes and of the maximum of the Aao correction. Within the 
presently allowed range of a, the largest difference in |Aao | between the exact 
Aao correction from eq.(|I~7D and the first order approximation of eq.(^) shows 
up ( see fig. 0(c)) at a ~ 60° and a ~ 105° approximately, while it is negligible 
for a values close to tt/2. It is worth noticing that we have used a positive value 
for P/T as provided by the model. When reversing its sign, the correction Aa 
changes sign too. 

Rewriting the (second order) expressions for the two parameters ao and &o in 
terms of a , S and A P /A T in eqs. (|T7D and ( |T8| ) as : 

a + sin 2a — 2 (A P /A T ) cos 5 ( a cos a + sin a ) + a (A P /A T ) 2 = (23) 
b - 2 (A P /A T ) ( b cos 5 cos a - sin 5 sin a ) + b (A P /A T ) 2 = (24) 



and inserting A P /A T from ( p0|) into eqs. (j23|) , (gj) with the value of P/T provided 
by the model, we get two equations in the two unknowns a and 5, both a function 
of the (3 angle : a = a(/3) and 8 = S(/3) . 
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Figure 2: (a) current experimental limits in the (p, 77) plane ; (b) correlation 
between sin 2/3 and sin 2a ; (c) current bounds on A P /A T plotted vs. a for 
P/T = 0.055 . 
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Figure 3: Aa correction plotted against the weak phase a for 5 = and 
four different values of j3 — 11° (dashed upper curve), 18° (solid), 25° (dashed), 
32° (dotted) (a) : exact correction ; (b) : approximate correction to first order 
in A P /A T ; (c) : difference (a)-(b). 



10 



Once (3 is known from an independent measurement ( e.g.: from the decay 

(-) 
B° d — ► J/ipKg ) and the two parameters clq and &o are fitted from the B°— > 7i + n 

time-dependent decay rates, we can in principle extract the correct value of a, 

from the data, together with a measurement of the relative strong phase S. 



4 Numerical solutions for Aa and 5 

For a given value of /3, we solve numerically the system of two equations (|23| ) and 
(p4[) in the two unknowns a and 5 and find the correction Aa to be applied to 
the measured value a M of a defined as : 

Aa = a M — a (25) 

where a M is inferred from the fitted value of the a® coefficient, with no penguin 
corrections, via the relation : 

— sin 2a M = a (26) 

Since the value of the measured asymmetry parameter a® is related to the weak 
phase a M via the above circular relation, one cannot distinguish between the case 
where a M falls into the [45°, 135°] interval or outside. Therefore, for — 1 < a < 1 , 
two numerical solutions for the true value of a are found. The first one corre- 
sponds to the above interval for a M , while the second solution belongs to the 
[ 0°, 45°] interval for a < and to the [ 135°, 180° ] interval for a > . 

Before taking into account the constraints from the present experimental 
bounds on a and /3, we first examine the predictions for both solutions i.e.: 
with a spanning the full [0 , n — f3] interval. 



The penguin correction Aa 

The "penguin- corrected" weak phase a and the correction Aa resulting from 
our analysis are plotted in fig. £|(a) and fig. f|(b), respectively, against the mea- 
sured ( uncorrected ) value a M in B° — > ir + 7i~ for four different values of j3 — 
11°, 18°, 25°, 32° and with a fixed value of b = -5.0 1(T 3 . This value has been 
chosen as an example where the value of the bo parameter is very small, but not 
identically zero, to exemplify a case where 5^0. 

From inspection of fig. |](a) and fig. |](b) we note that : 

• the Aa correction increases when j3 decreases. A maximum Aa correction 
of about 16.5° is found for j3 = 11°. The correction does not exceed 10° for 
(3 values greater than 18°, approximately. 
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Figure 4: For four different values of (3 = 11° (dashed), 18° (solid), 25° (dashed), 
32° (dotted) and P/T = 0.055 : (a) "a corrected" vs. "a measured" ; (b) Aa 
correction vs. "a measured" ; (c) relative strong phase 5 vs. "a measured" ( see 
section 4 for a detailed explanation of this picture ). 
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• maximal Act corrections for different j3 values occur when a M lies within 
the approximate interval ([70°, 100°]) and the true value of a corresponding 
to the maximum of the correction is found to increase when (3 decreases. 
However, to establish the maximal correction - allowed within the present 
limits on a and (3 - one should take into account the correlation between 
the allowed ranges of the two weak phases, as already pointed out for fig. [L| 

The "penguin-corrected" value of the asymmetry —sin 2a in B° d —* 7r + 7r~, corre- 
sponding to a M in the [45°, 135°] interval, is plotted in fig. || against the measured 
value of ao for four different values of (3 = 11°, 18°, 25°, 32° and with bo = 0. As 
expected from fig. |3](a) , we note that the measured value of a is significantly 
higher than the "true" value — sin 2a and that the correction increases for de- 
creasing (3 values. On the contrary, when a M lies outside the above interval, the 
correction ( not shown in this picture) to the true value of a is marginal and of 
opposite sign with respect to the previous case ( i.e.: the measured value of a^ 
is slightly larger than the "true" value ). 
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Figure 5: — sin 2a vs. measured a for four different values of (3 = 11° (dashed), 
18° (solid), 25° (dashed), 32°(dotted) with b = and P/T = 0.055 . 
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The strong phases difference 5 

In the above example, the values of the solution for the strong phases difference 
5 corresponding to a measured value a M are plotted as the family of curves in 
fig. |](c) for the same values of j3 used in the previous plots. In the picture, we 
note that the relative strong phase S shows a broad minimum which covers most 
of the region where a is bounded at present. In our example, the value of the 
b parameter being very close to zero, the minimum value of 5 turns out to be 
quite small (< 2° , approximately). On the contrary, the value of 5 shows a fast 
increase when a approaches the two "geometrical" limits or {n — f3). We want to 
stress that the dependence of 5 on a M in fig. f|(c) has no direct physical meaning. 
Instead, it describes the parametric behaviour of the solutions of (p3|) and ([24]) 
for the unknowns 5 and a when the value of the bo parameter is kept fixed at a 
non zero value and a is varied. The behaviour of 5 on the ascissa in fig. |](c) 
can be understood by considering the following (exact) relation : 

.* b [ C os2a-(A P /A T ) 2 ] 

tano = — : — — . t ,», (27) 

a [l + (A P /A T ) 2 }+sm2a v ; 



which can be easily derived from eqs.(|l7D and flTq) . 

The ratio A P /A T is found, from eg. (|20"D , to vanish in the limit a — > (it — (3), while 

A P /A T = P/T = 0.055 for a = (see also fig. |l|). Correspondingly, Aa — > for 

both a — » and a —>■ (n — 0). Therefore, when a is close to zero, the numerator 

of eq.(^7|) is almost constant, while the denominator approaches zero since, in 

this case, the Aao correction is almost negligible and the measured asymmetry 

oq does not differ significantly from the true asymmetry which is identically zero 

(see fig. |(b)). 

A similar, but not identical behaviour, takes place for a — > (7r — /3). This accounts 

for the fast increase of 5 on the final portion of the curves in fig. ||(c), drawn for 

the four different values of /3 when a approaches the respective values n — ft. 



4.1 The role of the 6q parameter 

The relation between 6 an d S, for given values of a and /3, is given in eq.(|l! 
Using the value of P/T from section 2, the numerator of fll8D is found to dominate 
over the weak cos 5 dependence of the denominator and therefore bo turns out to 
be approximately proportional to — sin 5 . This dependence is shown in fig. |6](a) 
where, taking (3 = 18° and letting — n < 5 < n, we plot bo as a function of 5 for 
two values of a = 104.8° (solid line) and a = 56.5° (dashed). In the same figure, 
the values of S corresponding to a minimum or a maximum for b are not too 
far from ±7r/2 for both curves in our example, as the absolute value of the bo 
parameter is maximal for cos<5 = 2(Ap/At) cos a / [1 + (Ap / At) 2 }- 
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Figure 6: With f3 = 18° in the four pictures : (a) b vs. S and (b) Aa vs. 
b for a = 104.8° (solid) and a = 56.5° (dashed). Taking four different values 
of b = -0.075 (dashed), -0.050 (solid), +0.025 (dashed), +0.075 (dotted) : 
(c) correction Aa Q vs. the measured asymmetry a Q ; (d) relative strong phase 5 
vs. "a measured". 15 



In order to show the dependence of the "penguin-correction" Aao , defined in 
eq . (|2l|) , on the actual value of the &o parameter, we stick to the above example 
and plot Aa as a function of 60 i n fig- H|(b). For — n < 5 < n , we found that 
Aao and bo satisfy an elliptic constraint (whose geometrical parameters depend 
on the choice of a and (3) . This is readily verified using the first order expansion 



( |22|) in A P /A T of equation (|T7|) and the analogous expansion of (|T8|). In this 



approximation, the constraint is simply 



r ^2 i2 , r b 12 _ l , 2i 

[ 2 (A P /A T ) sin a cos 2a l [ 2 (A P /A T ) sina 1 ' { 



Using instead the (exact) equations ( |17D and (|18D , the constraint is no longer 
represented by the approximate implicit form (ffijj ) and the ellipse is shifted along 
the Aa axis as in fig. ||(b) by a quantity which depends both on a and (3 (see 
next paragraph). 

The position of a point along the ellipse is parametrized in terms of the relative 
strong phase 5. For both curves in the example of fig. |6](b), a maximal positive 
(negative) correction Aao is reached for 5 = (5 = ±7r ) where 60 — 0. For not 
too large values of 8, which correspond to |6 I values close to zero, the correction 
Aao is positive and almost independent of 1 60 I - O n the contrary, the asymmetry 
correction Aa decreases rapidly as a function of \b | ( see fig. ||(b) ) when this 
parameter approaches its maximum allowed value (for a given a and (3) and 
the corresponding 5 is close to ±7r/2. For even larger strong phases, Aa may 
become negative. However, as pointed out previously, large values of \b \ would 
only occur for unexpectedly large strong phases. 

An example where the magnitude of 60 is kept small is shown in fig. ^|(c) where 
again we keep (3 = 18° fixed and, using the solutions of (^3|) and (0), we plot 
Aao vs. a for four different values of 60 =—0.075, —0.050, 0.025, 0.075 . 
In this range, we find that Aao scales approximately as cos S and therefore Aao 
is only marginally affected by 60 • I n this example, the corresponding solutions 
for 5 are those of fig. ^|(d) where 5 is plotted against the measured angle a M . 
In conclusion, the value of the 60 parameter is found to provide useful information 
on the strong phases difference 5. If 5 turns out to assume large values ( contrary 
to the conventional expectation ), then the parameter 6 is expected to be large 
in magnitude and to contribute significantly to the assessment of the "penguin 
corrected" value of a, which is otherwise determined by the value of a only. 
The current experimentally allowed range for the b Q parameter is discussed in 
the following paragraph. 
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Figure 7: Current bounds from CKM limits on : (a) maximal allowed value 
of the b parameter plotted against a ; (b) penguin-corrected asymmetry vs. 
measured asymmetry a ; (c) strong phase 5 vs. measured asymmetry a for 
b = ±0.1. (d) same as (c) for 6 = ±0.05 . 
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4.2 Bounds from the present CKM limits 

The maximum value |6 maa; | allowed by the present experimental CKM limits 
HI ! HZ! ) HH is plotted in fig. |7|(a) as a function of the true value of a when both 



a and (5 are allowed to vary through their common dependence on p and 77 within 
the domain shown in fig. 0(b). 

The values allowed at present by the experimental CKM limits for the "penguin- 
corrected" asymmetry — sin 2a are plotted in fig. 0(b) against the measured value 
of the asymmetry ao = —2 sin a M by varying — tc < 5 < it in such a way to 
account for all possible allowed values of b . As pointed out previously, the mea- 
sured value of ao is significantly higher than the "true" value of the asymmetry 
when a M belongs to the interval [45°, 135° ]. On the contrary, we can see in this 
picture that a M solutions located in the approximate range [20°, 45°] produce 
a marginal correction to the true value of a of opposite sign with respect to the 
previous case (i.e.: the measured value of a is slightly smaller than the "true" 
value ). 

Following a different approach, if we keep the magnitude of \b \ fixed and let 
a and (3 vary within the present CKM bounds, we can extract, as a solution 
of (p3|) and (|24]), the value of 5 which corresponds to a given observed value of 
the asymmetry ao . This is plotted in fig. 0(c), where we have allowed for our 
ignorance of the sign of b , and we have superimposed on the same picture two 
examples. The bounded area with positive (negative) values of 5 corresponds to 
&o = —0.1 and bo = 0.1, respectively. For large negative values of ao on the left 
side of figs. 0(c), (d) two solutions for 5 are found : the one with larger absolute 
values of 5 again corresponds to the approximate interval [ 20°, 45° ] for a M not 
(yet) ruled out by the present experimental limits. A similar example is shown 
in fig. |^(d) for 6 — ±0.05 and where \5\ is of order 10° over a large fraction of 
the accessible range of a . 

Next, we want to establish the combined bounds for the ao and b parameters. 
For a given choice of a and j3, the two parameters ao and bo are constrained 
onto one ellipse in the (ao , &o ) plane and the value of the parameter 5 is used to 
define the position of the point along the curve. From the (exact) equations (O) 
and (li"8|) , the constraint turns out to be an ellipse displaced along the ao axis by 
the amount ao : 

C-^f + (^L-? = 1 (29) 

where ao = [(A P /A T ) 2 — 1] sin la / C and the minor semi-axis is given by the 

expression : m = 2 (A P / A T ) sin a[(A P /A T ) 2 — cos 2a)] / C , while the major 

semi-axis |6 ma:I '| is given by : ( Jj Q max ^ 2 = 4(A P /A T ) 2 sm 2 a / C and the common 

denominator is : C ~ 1 — 2(A P /y4 T ) 2 cos2a. 

By varying a and (3 within their allowed domain in the (p, rj) plane, a closed 

boundary in the (a , bo ) plane is obtained which is shown in fig. [| 

This boundary can also be seen as the projection onto the (ao , &o ) plane of the 



surface represented in fig. [| In this picture, we plot the value of sin 5 extracted 
as a solution of ( p3|) and fl24j) for each point in the allowed region of the (a , &o ) 
plane with the assumption that the corresponding angle /3, as well as the other 
two measured quantities a® and 60 , are perfectly known, i.e.: no experimental 
uncertainties are introduced at this level. 

As a final remark, we want to underline that the method used to derive the bounds 
presented in this section, although applicable in general, depends numerically on 
the value of P/T = 0.055 from section 2 and, in this respect, the actual values of 
the bounds have to be considered as "model-dependent" . 
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Figure 8: Combined bounds for a^ and b from the present CKM limits and 
with P/T = 0.055 . 
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5 Conclusions 

We have described a quantitative procedure to extract ( in a model-dependent 
way ) the value of the CKM angle a from a measurement of B° — > 7r + 7r~ in the 
presence of penguin pollution. We have shown that uncertainties of strong inter- 
action origin can be controlled and that a model-dependent penguin correction 
to the measured value of a from B° — »■ 7r + 7r~ can be evaluated and parameterized 
in terms of the CKM angle (3 and of the parameters a$ and &o extracted from a 

time-dependent analysis. Also, the cos Ami term in a time- dependent analysis 

(-) 
of B°— > 7T + 7r is found not to be sufficient to provide a direct measurement of 

the size of the penguin, but to provide instead useful information on the phases 

of strong origin. 

The model-dependence of the procedure stems from the numerical dependence 

on the ratio of penguin-to-tree hadronic matrix elements P/T (see section 2) of 

the penguin correction Aa (to the weak phase a), of the corresponding correction 

b 
Aa (to the observed asymmetry a ) and of the ratio -A*. 

The model-dependent procedure to extract a and 5 proposed in this paper is 
applied to one final numerical example where we take the value P/T = 0.055 
from section 2 and assume for (3 a measured central value of (3 = 18°. An 
observed asymmetry a$ ps —0.14 would, in this example, require a large pen- 
guin correction Aao ~ 0.32 and the extracted "true" asymmetry value would be 
ao corrected Ki —0.46. In terms of a, a correction Aa pa 10° should be applied to the 
measured central value a M ~ 86° and would result into a "penguin corrected" 
value of a ~ 76°. An obvious comment is that a Aa correction of order 10% 
corresponds to a significantly larger correction in terms of Aa as a consequence 
of the — sin 2a dependence of the asymmetry on a. In our example, a measured 
central value of bo = —0.1 would correspond to a value for the relative strong 
phase 5 ~ 16°. 

The experimental errors on the parameters ao and bo from a time-dependent 
analysis of B° — ► 7r + 7r~, together with the expected error from an independent 
measurement of sin 2/3, propagate into an uncertainty on the asymmetry correc- 
tion Aa and on the extracted value of S. 

The evaluation of such uncertainties as well as the assessment of the experimental 
sensitivities to a and 5 is beyond the scope of the present paper. 



Acknowledgments 

We would like to thank P.Colangelo for useful discussions. 



21 



References 

[I] For a review see, e.g., Y. Nir and H.R. Quinn, Annu. Rev. Nucl. Part. Sci. 
42 (1992) 211, and references therein. 

[2] M. Gronau, Phys. Lett. B300 (1993) 163. 

[3] M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34 (1987) 103. 

[4] L. Wolfenstein, Phys. Rev. D 43 (1991) 152; and references there. 

[5] G. Kramer, W.F. Palmer and Y.L. Wu, Report DESY 95-246 (1995). 

[6] H. Simma and D. Wyler, Phys. Lett. B272 (1991) 375. 

[7] A.J. Buras and R. Fleischer, Phys. Lett. B341 (1995) 379. 

[8] M. Bander, D. Silverman and A. Soni, Phys. Rev. Lett. 43 (1979) 242. 

[9] J.F. Donoghue, E. Golowich, A. Petrov and J.M. Soares, Phys. Rev. Lett. 77 
(1996) 2178. 

[10] A.N. Kamal, Int. J. Mod. Phys. A7 (1992) 3515. 

[II] B. Blok and I. Halperin, Phys. Lett. B385 (1996) 324. 

[12] M. Gronau and D. London, Phys. Rev. Lett. 65 (1990) 3381. 

[13] C. Hamzaoui and Z. Xing, Phys. Lett. B360 (1995) 131. 

[14] N.G. Deshpande and X.-G. He, Phys. Lett. B384 (1996) 283. 

[15] G. Kramer and W.F. Palmer, Phys. Rev. D 52 (1995) 6411. 

[16] A. Deandrea, N. Di Bartolomeo, R. Gatto, F. Feruglio and G. Nardulli, 
Phys. Lett B320 (1994) 170. 

[17] J. P. Silva and L. Wolfenstein, Phys. Rev. D 49 (1994) 1151. 

[18] Z. Xing, Nuovo Cim. A108 (1995) 1069. 

[19] N.G. Deshpande and X.-G. He, Phys. Rev. Lett. 75 (1995) 1703. 

[20] M. Gronau, O.F. Hernandez, D. London and J.L. Rosner, Phys. Rev. D 52 
(1995) 6356. 

[21] A.S. Dighe, M. Gronau and J.L. Rosner, Phys. Rev. D 54 (1996) 3309. 

[22] R. Fleischer and T. Mannel, Report TTP 96-49 (1996). 



22 



[23 
[24 
[25 

[26 

[27; 
[28 
[29 

[30 

[31 
[32 

[33 

[34 
[35 
[36 

[37; 

[3 



A.J. Buras and R. Fleischer, Phys. Lett. B360 (1995) 138. 

F. DeJongh and P. Sphicas, Phys. Rev. D 53 (1996) 4930. 

A.E. Snyder and H.R. Quinn, Phys. Rev. D 48 (1993) 2139; H.J. Lipkin, Y. 
Nir, H.R. Quinn and A. Snyder, Phys. Rev. D 44 (1991) 1454. 

R. Aleksan, F. Buccella, A. Le Youanc, L. Oliver, O. Pene and J.-C. Raynal, 
Phys. Lett. B356 (1995) 95. 

H. Fritzsch, D. Wu and Z. Xing, Phys. Lett. B328 (1994) 477. 

D.Du and Z. Xing, Phys. Lett. B280 (1991) 292, and references therein. 

A.J. Buras, M. Jamin, M.E. Lautenbacher and P.H. Weisz, Nucl. Phys. B370 
(1992) 69; ibid. B400 (1993) 37; A.J. Buras, M. Jamin and M.E. Lauten- 
bacher, Nucl. Phys. B400 (1993) 75. 

M. Ciuchini, E. Franco, G. Martinelli and L. Reina, Nucl. Phys. B415 (1994) 
403. 

N.G. Deshpande and X.-G. He, Phys. Rev. Lett. 74 (1995) 26. 

M. Gronau, O.F. Hernandez, D. London and J.L. Rosner, Phys. Rev. D 52 
(1995) 6374. 

M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34 (1987) 103. 

Review of Particle Properties, Phys. Rev. D 54 (1996) 1. 

Y. Koide, Shizuoka report US-94-05 (1994). 

A. Ali and D.London, Report DESY 96-140 (1996). 

A. Ali, Report DESY 96-248 (1996). 

J. Rosner, Report EFI-96-46 (1996). 



23 



